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Self-diffusiophoresis of Janus particles in near-critical mixtures 


Alois Wiirger 

Laboratoire Ondes et Matiere d’Aquitaine, Universite de Bordeaux & CNRS, 

351 eours de la Liberation, 33405 Talenee, Franee 

We theoretically study the self-propulsion of a laser-heated Janus particle in a near critical water- 
lutidine mixture, and relate its velocity and squirmer parameter fd to the wetting properties of 
its two hemispheres. For non-ionic surface forces, the particle moves the active cap at the front, 
whereas a charged hydrophilic cap leads to backward motion, in agreement with experiment. Both 
Vp and /3 show non-monotonic dependencies on the heating power, and may even change sign. The 
variation of /3 is expected to strongly affect the collective behavior of dense squirmer systems. 
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In recent years, artificial microswimmers have been 
realized by Janus particles which move along the con¬ 
centration or temperature gradients generated by their 
own chemical or thermal activity Oriented au¬ 

tonomous motion has been achieved throng dynamical 
feedback [ 6 ] or rectification in a periodically structured 
channel [7], opening applications such as targeted trans¬ 
port and pumping of passive particles. In dense systems, 
active Janus particles aggregate in dynamical clusters 
019]. This observation was related to short-range hydro- 
dynamic effects HDHH], in terms of the squirmer model 
originally developed for the motility of bacteria. Self¬ 
propulsion mechanisms have generally a strong diffusio- 
phoretic component nmn]; in the case of ionic molecular 
solutes, self-generated electric fields and ion effects may 
contribute to the motion [laiis]. 

Diffusiophoresis was first rationalized by Derjaguin et 
al. [17], when observing that wax particles dispersed 
in a non-uniform glucose solution, migrate toward lower 
sugar concentration. Because of its unfavorable interac¬ 
tion with wax > 0 ), sugar is depleted in the boundary 
layer, the adsorption parameter F = — 1) 

is negative, and the wax particle migrates to lower glu¬ 
cose content at the velocity [18] 

Vp = -^FVn, (1) 

where 77 is the solvent viscosity. For a molecular solute 
that is attracted by the surface, one has F > 0 and the 
particle moves toward higher concentration. 

There is no such simple rule for self-propelling Janus 
particles, where both the adsorption parameter and the 
concentration gradient vary along the particle surface, 
and where catalytic coating may result in a multicom¬ 
ponent boundary layer na ES]. A particularly intri¬ 
cate situation occurs for hot Janus particles in a near- 
critical water-lutidine mixture [4], which migrate in a 
self-generated composition gradient. Though their mo¬ 
tion is clearly related to the wetting properties of their 
active and passive surfaces 0120 ], there is at present no 
explanation for the sign and magnitude of the velocity. 

In this Letter we study self-diffusiophoresis in such 
near-critical binary mixtures 0 0 HO]. Starting from 




FIG. 1: Critical droplet (0 > (f)^) surrounding a particle with 
hydrophilic surface coating on both the cap (red) and the un¬ 
capped hemisphere (white). The surface temperature profile 
T(r) (calculated from the thin-cap limit of Ref. [21]) and the 
composition (from with C = 100 °C) are plotted as a func¬ 
tion of the cosine of the polar angle. The scale of the reduced 
temperature (T — To)/(Tc' — To) is on the left ordinate, and 
that of the water content 4> — 4>c on the right one. The pa¬ 
rameter T = (Tm — To)/(Tc' — To) depends on the temperature 
at midplane Tm- 


the properties of the demixing volume surrounding the 
particle, we evaluate the velocity Vp and the squirmer pa¬ 
rameter {3, both of which depend in an intricate manner 
on the heating power and the dispersion forces exerted 
by the two hemispheres. Finally, we compare with recent 
experiments and discuss charge effects. 

The critical droplet. Fig. 1 illustrates an active Janus 
particle in a water-lutidine mixture at the critical wa¬ 
ter content (I)q = 0.72 and at a bulk temperature Tq 
which is slightly below the critical value Tq = 34.1 °C. 
Illuminating the particle with a laser beam, results in 
a temperature profile T(r) that exceeds Tc on part or 
all of the surface, and thus causes local demixing. As¬ 
suming a quadratic relation to the local composition, 
T — Tc = C{(l) — (l)cY^ one finds the change of water 
content 


= ( 2 ) 

where the two signs correspond to water-rich and 
lutidine-rich phases, and where C ^ 100 °C [ 22 ]. This 
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FIG. 2: Demixing volume surrounding a particle with hy¬ 
drophobic cap and hydrophilic remainder. Heating power in¬ 
creases from a) to c); the right panels show temperature and 
composition prohles as in Fig. 1. a) The critical condition is 
satisfied only in a thin lutidine-riche droplet that partly cov¬ 
ers the cap. b) Most of the surface is above Tc; the cap is 
surrounded by the lutidince-rich phase, and the hydrophilic 
remainder by the water-rich phase. The phase boundary is 
modelled by a factor tanh(cos ^/co) with co = 0.1. c) Strong 
heating leads to T > Tc on the entire surface and to a spher¬ 
ical demixing volume. 


mean-field relation ceases to be valid at the critical point 
where composition fluctuations become long-range. 

A more complex situation occurs for a particle with 
a hydrophobic cap, as illustrated in Fig. 2. If Tc is 
reached on both hemispheres, the critical volume splits 
in lutidine-rich and water-rich compartments, with the 
phase boundary attached at the particle’s midplane [23] . 
This separation is possible since the interface tension be¬ 
tween the spinodal phases, 7 < 10“^ N/m for T < 36 °C 
[ 22 ] . is much smaller than the particle’s surface energy. 
For the strong-heating case (c), the critical volume is al¬ 
most spherical, and its radius ra much larger than that 
of the particle, a. Experiments cover the whole range 
from a small demixing area to a large critical droplet of 
tens of microns [ 20 ] . 

Slip velocity. Contrary to the demixing volume, the 
boundary layer is not in a quiescent state, but shows a 
non-uniform pressure and a steady diffusion current. As¬ 
suming that mutual diffusion of water-lutidine is faster 
than advection and inserting the current in Stokes’ equa¬ 


tion, one finds the effective slip velocity [24] . 


ksT ^dcj) 
VT] dx 


(3) 


where we have defined the mean inverse molecular volume 
v~^ = (j)cvY^-\-{l — (j)c)v~^ and the adsorption parameter 


— p-^i 

r= / dzz- -^------. (4) 

The effective potential of water and lutidine, and 7 /^^ 
are given in units of the thermal energy, and vanish well 
beyond the interaction range A. In the dilute limit 0^0 
and with u = ~ 1)1 recovers Derjaguin’s 

adsorption factor in (Ih. Note that that the slip velocity 
Vs does not depend on the parallel force component 
but on the composition gradient only [28] . 

For an order-of-magnitude estimate, it is convenient to 
explicit the adsorption parameter for a square well poten¬ 
tial of width A and prefactor -0, where the integral in Q 
gives ^A^ with -0^ ^ in the second factor. A strongly 
hydrophobic surface repels water and attracts lutidine 
('0^ > such that F < 0, whereas a hydrophilic sur¬ 
face is characterized by F > 0. With typical parameters 
m we find F ^ 10 m^ and, supposing an interaction 
length A of a few A, we deduce r\j 10 Then the ad¬ 
sorption parameter F = takes a constant 

value Fcap on the cap, and Func on the remainder. 

Self-propulsion. The particle velocity is obtained by 
averaging the slip velocity over the surface, = — (v^) 
m- For an axisymmetric particle one finds 

Vp = ^ 0 ^ f dc{l - c^)Tdc<p, (5) 

2vr]a 

where we have used the relation dx = a~^y/l — c^dc be¬ 
tween the local coordinate x and the cosine of the polar 
angle c = cos 0. Sign and magnitude of the velocity are 
determined by the product of the adsorption factor F and 
the derivative of the composition gradient dcf). A par¬ 
ticularly complex behavior occurs for cases as in Fig. 2, 
where Fcap and Func take opposite signs and where the 
gradient dc<f is largest in the midplane area. 

In Fig. 3 we plot the particle velocity as a function 
of the surface temperature at midplane T^, in terms of 
the reduced quantity r = (T^ — Tq)/{Tc — Tq). Criti¬ 
cal conditions on the summit of the cap are reached at 
Tc = oa the particle’s midplane at = 1, and 

on the entire surface at = 1 + 1/v^- The behavior 
of the velocity is to a large extent determnined by the 
ratio of adsorption parameters f = Tunc/Tcap; that is, 
by the wetting properties of the two hemispheres. We 
distinguish three parameter ranges. 

(i) For Tc < T < 1, the critical droplet covers only 
part of the cap, as in Fig. 2a. Since in this range, F 
and dc(f carry the same sign, the particle moves forward 
{vp > 0) for both hydrophobic and hydrophilic coating. 
The velocity is independent of Func- 
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FIG. 3: Self-propulsion velocity up as a function of r = (Tm — 
To)/(Tc' — Fo), for different values of absorption parameters. 
Positive Vp means that the particle moves forward (cap at the 
front). We fix the adsorption parameter of the cap and vary 
the ratio ^ = Func/Fcap- Self-propulsion sets in at re = 

If the critical droplet extends to both hemispheres (r > 1), the 
velocity strongly depends on ^ and changes sign for ^ < — 1. 
The velocity scale is of the order of /xm/s. 


(ii) In the range 1 < r and ^ > 0, both hemispheres 

contribute to the integral in The cusp at oc¬ 
curs because of the large derivative 9c0 on the uncapped 
hemisphere; see the right panel of Fig. 1. The veloc¬ 
ity goes through a maximum at r ^ 1.5, where most of 
the particle is covered by the critical volume; at strong 
driving r ^ 1, it increases as Up oc (1 + El]- 

(iii) For 1 < r and ^ < 0, the velocity is to a large ex¬ 
tent determined by the change of sign of (/) — at mid¬ 
plane, where the contributions of the two hemispheres 
partly cancel in ([^ . If the adsorption parameter is larger 
on the uncapped part, ^ < —1, it dominates the velocity 
and finally results in a change of sign; well beyond rs 
one finds Vp ex {1 — ^)\/r [24] . 

Several features can be traced back to the relation be¬ 
tween the laser intensity I and the excess temperature. 
With Tm — To = /xa/2/^, the heat conductivity and 
the absorption coefficient per unit area x, one finds 

Vp oc-^-, Ic =- [Tc - To). (6) 

In a very narrow range above rc one has f{x) = (in¬ 
visible in Fig. 3); the cusp at Tm follows the law yjI — Im- 
We briefly discuss the above result in view of recent ex¬ 
periments. At small or moderate driving we expect the 
particles to move the cap at the front. This agrees with 
observations on carbon-capped silica beads [9] and gold- 
capped particles with hydrophobic coating [20] . In addi¬ 
tion, the size dependencies of and the overall shape of 
Vp agree well with the data of Buttinoni et al. (Fig. 4a of 



FIG. 4: Squirmer parameter /3. For t ^ tc, the active area 
reduces to a small spot, resulting in = 5 [24]. The diver¬ 
gency for ^ = —3 occurs where the velocity Vp is zero. 

m), measured for beads of different radius (a = 0.5/im 
and 2.13/im) in the range r < 1 [24]. A strong discrep¬ 
ancy occurs for particles with hydrophilic caps, which 
move the cap at the rear, at both weak [4] and strong 
driving [20] . 

Some aspects of Fig. 3 disagree with a very recent the¬ 
ory paper m- As two main differences we note that Ref. 
[30] (i) assumes an isothermal cap with zero slip veloc¬ 
ity and (ii) considers quite large velocities Vp > lO/xm/s, 
where the water-lutidine kinetics is governed by advec- 
tion rather than diffusion [24]. 

Squirmer parameter The interaction of a microswim¬ 
mer with a wall and collective effects are to a large extent 
determined by the squirmer parameter (3 [TOUT^ . which 
is defined through the even component of the slip veloc¬ 
ity Vs = v^^ sin^(l + PcosO) [31]. A “puller” is propelled 
by the activity of its front hemisphere {P > 0), and a 
“pusher” by its back part (/3 < 0). 

In Fig. 4 we plot /3 as a function of r for differ¬ 
ent adsorption parameter ratios At the onset of self¬ 
propulsion, where the active area is reduced to a small 
spot at the summit of the cap, one finds 

/3 = 5 {j^Tc)- (7) 

(For a particle moving the active spot at the back, one 
has (3 = —5.) With increasing driving, (3 decreases rapidly 
and strongly depends on the reduced temperature r and 
the parameter Opposite affinities of the two hemi¬ 
spheres, C < 0, may result in pullers or pushers of vari¬ 
able strength; the singularity for ^ = —3 occurs where 
the particle velocity changes sign. This means that a 
tiny change in the driving could significantly modify p 
and thus the collective behavior [TQUT2] . 

Charge effeets. Surface charges have been shown to be 
relevant for the reversible aggregation of polystyrene par- 
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hydrophilic 



FIG. 5: Schematic view of the electric double layer with 
screening length . The ions within the demixing volume of 
thickness a diffuse toward higher water content. As a result, 
the particle moves backward (cap at the rear) for hydrophilic 
coating, and forward (cap at the front) for hydrophobic coat¬ 
ing. 


tides in a near-critical water-Intidine mixture [32], and 
they may even change the sign of the crtical Casimir ef¬ 
fect [33]. Here we complete the above discussion of self- 
diffusiophoresis by including charge effects, as a possible 
explanation for the backward motion of beads with hy¬ 
drophilic cap. 

A composition gradient along a charged surface 
gives rise to two distinct effects: the drift of the mobile 
counterions due to the ion-specific thermodynamic force 
—V/i, and the non-uniform properties of the electric- 
double layer, very much like the thermal forces in a tem¬ 
perature gradient [28] . Here we discuss the ion-drift term 
only, and reduce the chemical potential to the electro¬ 
static self-energy of a monovalent ion of radius a^, 


yu/ — . i u j 

Sneam 

The variation of the composition (j) in the demixing vol¬ 
ume is illustrated in Fig. 5. The dependence of the 
permittivity s on (p gives rise to a thermodynamic force 
density —pdxli with the ion concentration p. This ion 
current drags the fluid along the particle surface and thus 
induces a slip velocity 

1 

Vs = - - / dzzpdxp. (9) 

V Jo 

Spelling out the gradient, dxP = —p{d(f)\ns)dx(l), assum¬ 
ing the linear law 5 = + (1 — <P)£i for the permit¬ 

tivity of water-lutidine, and using ^ one finds 

The particle velocity Vp is given by the surface aver¬ 
age of the negative slip velocity. Evaluating the counte¬ 
rion concentration in Debye-Hiickel approximation, p = 
(5/e)|C|e“^^, we find [24] 


= __£lC|_ P (1 - c^)dc4> 

Sirriama Ti [1 + (cr/t)-!]^ ’ 


( 10 ) 


which strongly depends on the ratio of the screening 
length and the thickness a of the demixing volume. 


With typical parameters the prefactor takes a value of 
millimeters per second. Taking cr/i: ^ ^ and dc<p ^0.1 
as suggested by Figs. 1 and 2, we find Vp ^ /im/s, which 
corresponds to measured values. So far we have consid¬ 
ered the salt-free case where ^ hundreds of nanome¬ 
ters. Adding salt would result in diffusiophoresis in a 
non-uniform electrolyte j34f[36] . Moreover, we have dis¬ 
carded specific-ion effects which are not small in general 

[HIST]. 


According to (10), charged particles move in the direc¬ 
tion opposite to the composition gradient, that is, cap at 
the rear for water-adsorbing (hydrophilic) coating. This 
is precisely what was observed in experiments on parti¬ 
cles with hydrophilic gold caps nnn]. The ionic end- 
groups used in these studies (11-mercapto-undecanoic- 
acid) cause a ^-potential of about —50 mV [38], which 
results in a negative velocity Vp of microns per second. 


On the other hand, the cap-at-the-front orientation, 
expected for lutidine-adsorbing (hydrophobic) charged 
surfaces, is probably of little relevance: The carbon caps 
[9] and hydrophobic gold caps gold caps (functionalized 
with 1-octadecanethiol) [20] carry only weak charges; as 
a consequence, their forward motion is due to the disper¬ 
sion forces underlying <§ 

Conclusion. Hot Janus particles in a near-critical 
water-lutidine mixture move due to their self-generated 
composition gradient. Our analysis reveals two main 
mechanisms: The dispersion forces exerted on water and 
Intidine result in a positive velocity Vp > 0^ whereas 
at a charged surface, the counterions migrate toward 
higher water content and thus drive hydrophilic parti¬ 
cles backward, Vp < 0. The first effect accounts for the 
observed forward motion of uncharged particles [9| |20] , 
and the second one for the backward motion of beads 
with charged hydrophilic caps [1[20]. 


Both the velocity Vp and the squirmer parameter (3 
show a non-monotonous variation with the wetting prop¬ 
erties and surface temperature; this implies that the hy¬ 
drodynamic interactions related to {3 strongly depend on 
the driving power. The phase behavior of squirmer sys¬ 
tems is very sensitive to the value of p [TOUT^ : in view 
of Fig. 4 it could be changed by tuning the heating. 


As an outlook, our findings, both on the wetting prop¬ 
erties and on charge effects, could be relevant for other 
driving mechanisms of active particles. Moreover, they 
suggest that the overlap of the critical droplets of nearby 
Janus particles should give rise to a complex interaction 
pattern, which could affect the observed aggregation be¬ 
havior and result in a variety of reversible ordered states, 
similar to those realized recently with Janus particles in 
a homogenous near-critical water-lutidine mixture [23] . 
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